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ABSTRACT 


The  simplest  version  of  the  Lauritzen-Hoffman  (LH)  model  of  polymer 
crystallization  which  applies  to  infinitely  long  model  polymer  molecules 
crystallizing  on  an  existing  substrate  of  infinite  width  is  reexamined.  The 
mathematical  expressions  for  the  model  free  energy  barriers  are  observed  to 
take  on  negative  values  at  high  supercooling.  Since  such  negative  barriers 
appear  to  be  physically  unrealizable  for  the  crystallization  process,  the  LH 
model  is  extended  only  by  imposing  a  mathematical  constraint  on  the 
expressions  for  the  barriers  to  forbid  them  from  ever  being  negative.  The 
extended  model  contains  one  parameter  7  which  varies  from  zero  to  one  and  is 
analogous  to  the  parameter  ^  of  the  LH  model.  For  all  values  of  7  less  than 
one,  the  extended  model  predicts  a  finite  lamellar  thickness  at  every 
supercooling;  moreover,  this  thickness,  at  large  undercooling  decreases 
monotonically  with  increasing  undercooling  in  agreement  with  experiment  but  in 
marked  contrast  to  the  LH  model  which  exhibits  the  well-known  S£  catastrophe. 
The  relative  insensitivity  of  the  calculated  lamellar  thicknesses  to  the 
parameter  7  supports  the  use  of  low  parameter  values  such  as  zero  as  a  first 
approximation  for  mathematical  convenience  in  practice . 
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I .  INTRODUCTION 

Recently,  the  isothermal  (unoriented)  crystallization  of  poly (vinylidene) 

fluoride  (PVF2)  from  the  melt  in  the  presence  of  a  high  static  electric  field 

has  been  studied  experimentally  and  theoretically  as  a  simple  model  system  for 

the  investigation  of  the  structure -property  relationships  of  this 
1112 

polymer.  •  ’  •  Of  the  four  well-known  crystalline  forms  a,  7,  and  S  (or 

II,  I,  III,  and  IV)  of  PVF2:  however,  the  phase  which  has  the  greatest 
potential  for  applications  via  its  large  spontaneous  polarization  '  is  not 
produced  on  crystallization  from  the  melt.  This  phase- -namely  ^--has  been 
observed  to  exhibit  highly  oriented  growth  (mixed  with  some  unoriented  a  phase 
growth)  during  the  initial  stages  of  crystallization  from 

solution  under  a  high  static  field  of  highly  plasticized  PVF2  (i.e.  of  a 

concentrated  solution  of  tricresyl  phosphate  and  PVF2  in  which  the  tricresyl 

phosphate  content  decreases  by  evaporation) ;  at  longer  times ,  the  formation  of 

I  3 

the  nonpolar  a  phase  predominates.  '  The  decrease  in  plasticizer  content  and 
the  subsequent  crystal  growth  behavior  implies  that  the  local  electric  field 
in  the  solution  region  changes.  As  part  of  the  continuing  effort  simply  to 
understand  the  structure -property  relationships  of  PVF2  and  because  of  its 
practical  importance,  our  ultimate  goal- -despite  the  complexity  of  the  system 
described- - is  to  develop  a  theory  or  model  which  can  account  for  its 
crystallization  behavior. 

As  in  the  case  of  isothermal  crystallization  of  a  and  7  phase  from  the 

I  2 

melt  in  an  electric  field,  ’  a  theory  of  isothermal  crystallization  of  a,  p, 
and  6  phase  from  concentrated  solution  in  an  electric  field  would  be  based  on 
"classical"  and  "polymer"  theories  of  nucleation  and  growth  in  the  absence  of 
an  applied  field.  Most  importantly,  the  nucleation  barrier  or  activation  free 
energy  barrier  for  nucleation  would  certainly  be  different  in  the  presence  of 
the  field  than  in  its  absence;  and  this  barrier  has  been  seen  to  be  of 
fundamental  importance  in  the  theories  of  polymer  crystallization,  the 
simplest  of  which  is  the  LM  or  Lauritzen-Hoffman  theory.  One  possibly 
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unrealistic  feature  which  seems  to  have  been  incorporated  into  this  theory  in 
order  to  simplify  it  is  that  the  nucleation  barrier  is  not  constrained  in  the 
theory  to  take  on  only  nonnegative  values.  The  word  "barrier"  connotes  a 
positive  quantity,  and  furthermore,  the  LH  theory  is  based  on  transition  state 
theory  in  which  the  barrier  corresponds  to  an  intermediate  configuration  or 
transition  state  of  the  system  which  is  at  a  free  energy  maximum  relative  to 
some  initial  and  final  state  of  the  system.^  Moreover,  the  LH  theory 
exhibits,  in  contrast  with  experiment,  the  Si  catastrophe  wherein  the 
calculated  average  lamellar  thickness  i  suddenly  passes  through  a  minimum  and 
becomes  infinite  at  a  temperature,  T^,  corresponding  to  a  moderately  large 
undercooling;  and,  in  fact,  the  nucleation  barrier  in  this  theory  is  positive 
for  all  T  >  T^,  is  zero  at  T  “  and  is  negative  for  all  T  <  for  the 

special  case  which  Lauritzen  and  Hoffman  have  recently  considered. 
Therefore,  prior  to  developing  an  extension  of  the  LH  theory  which  would 
involve  ascertaining  the  effect  of  an  electric  field  on  the  nucleation 
barrier,  we  try  to  extend  the  LH  theory  to  larger  undercooling  by 
incorporating  into  it  the  assumption  that  free  energy  barriers  cannot  be 
negative.  Note  that,  unlike  in  the  LH  theory  of  polymer  crystallization, 
barriers  in  classical  nucleation  theory  are  never  negative;  however,  the 
classical  theory  does  not  explicitly  take  into  account  polymer  chain  folding, 
and  for  that  reason,  we  have  not  yet  considered  modifying  the  Marand  and  Stein 
theory^ ■ ^  of  crystallization  from  the  melt  to  treat  the  PVF2/tricresyl 
phosphate  crystallizing  solution. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section  II,  the 
LH  model  is  described.  The  kinetic  treatment  of  the  LH  model  is  given  in 
Section  III.  The  rate  constants  needed  for  this  treatment  are  determined  in 
Section  IV.  Next,  our  extension  of  the  LH  model  is  described  in  Section  V; 
the  conditions  which  determine  the  sign  of  the  free  energy  of  formation 

of  that  portion  of  a  model  polymer  molecule  which  crystallizes  first  on  an 
existing  crystal,  are  found  in  Section  VI.  A  summary  of  the  expressions  for 
the  barriers  in  our  model  is  given  in  Section  VII  along  with  the  expressions 
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for  the  average  lamellar  thickness.  In  Section  VIII,  the  variable 
transformations  required  as  a  preliminary  to  numerical  integration  are 
introduced.  Results  and  discussion  appear  in  Section  IX,  and  conclusions  are 
given  in  Section  X. 


II.  THE  lAURITZEN-HOFFMAN  MODEL 

1  2 

The  model  to  be  extended  is  one  version  ’  of  the  well-known  Lauritzen- 

Hoffman  (LH)  model  of  polymer  crystallization.  Our  description  of  this 

version  is  as  follows.  The  model  polymer  molecules  are  assumed  to  be 

infinitely  long  and  crystallize  on  an  existing  crystalline  face  or  substrate 

which  is  assumed  to  be  infinitely  wide  (i.e.  the  fact  that  its  width  is  finite 

is  ignored)  .  A  sequence  of  length  f  of  polymer  segments  of  width  a  and 

thickness  b  as  well  as  the  volume  associated  with  that  sequence- -which  is 

taken  to  be  a  parallelepiped  of  length  i,  width  a,  and  thickness  b--is 

designated  as  a  stem.  Only  stems  of  length  i  can  cirystallize  on  an  existing 

face  of  length  f ,  but  the  length  i ,  the  lamellar  thickness ,  can  vary  from 

crystal  to  crystal.  (Of  course,  a  and  b  cannot  vary  from  face  to  face.)  Any 

sequence  of  length  i  of  segments  of  a  model  molecule  can  be  placed  first  on  a 

given  face  and  upon  placement,  is  designated  as  the  first  stem.  The  free 

energy  of  formation  of  the  first  stem  is 

-  A^q  ■  A^j^  -  0  or  A^j^  -  2aha'^  +  2bi<7  -  abiAf 

where  Af  >  0  is  the  free  energy  of  fusion  per  unit  volume  at  a  temperature  T 

below  the  melting  point  T^  of  a  crystal  of  very  large  2  (Af  -  0  at  T  -  T^) , 

where  a  is  the  lateral  surface  free  energy  per  unit  area  (i.e.  that  associated 

with  the  surfaces  of  area  bi  and  ai  of  a  stem) ,  and  where  is  the  surface 

free  energy  per  unit  area  associated  with  the  cilium  that  protrudes  through 

1-3 

each  of  the  surfaces  of  area  ab  of  the  first  stem.  Recently,  has  been 

assumed  to  be  zero;  generally,  one  can  have^  0  <  .  All  surface  free 

energies  per  unit  area  in  the  model  are  assumed  to  be  independent  of  T  and  2. 
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(See  Figure  2(a)  of  Reference  1  or  Figure  22  of  Reference  2.)  The  placement 
of  each  subsequent  stem  involves: 

1.  the  destruction  of  the  cilium  associated  with  one  of  the  surfaces  of 
area  ab  of  an  adjacent  stem  already  crystallized, 

2.  an  adjacent  reentry  and  the  formation  of  a  tight  fold  associated  with 
two  surfaces  of  area  ab,  and 

3.  the  formation  of  a  cilium  associated  with  the  remaining  surface  of 
area  ab  of  the  stem  being  placed. 

Only  adjacent  reentry  and  hence  only  tight  folding  is  incorporated  in  this 
version  of  the  model. 

The  free  energy  of  formation  of  the  vth  stem  (i/  >  1)  is  therefore 

^^‘^e  '  abiAf 
or 

A^  -  ,  -  2aba  -  abiAf  ■  -E 

where  A^^  is  the  free  energy  of  formation  of  a  group  of  u  stems  (relative  to 
A^q  ■  0)  and  where  is  the  surface  free  energy  per  unit  area  associated  with 

half  of  a  fold.  Iteration  of  A^^^  -  -  E  (i/  >  1)  gives 

A^^  -  A^j^  -  (v-l)E 

-  2bi<7  +  2aba'  -  2aba^  +  «/ab(2a  -iAf )  . 

6  0  6 

In  order  that  stem  additions  subsequent  to  the  placement  of  the  first  stem  be 
thermodynamically  favorable,  i.e.  in  order  that  they  would  in  fact  occur,  one 

must  impose  the  constraint  E  >  0  and  consequently  i  >  By  contrast,  A^^ 

can  be  positive,  zero,  or  negative;  E  >  0  guarantees  that  A^^  <  0  will  occur 
for  finite  i/.  Note  the  sign  conventions  for  Ai^j^  and  E. 
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III.  THE  KINETIC  TREATMENT  OF  THE  LADRITZEN- HOFFMAN  MODEL 

1  2 

Our  description  of  the  kinetic  treatment  *  of  the  LH  model  is  as 
follows.  The  following  assumptions  are  made: 

1.  Assume  that  transition  state  theory  can  be  utilized  to  describe  the 


kinetics  of  the  LH  model  of  polymer  crystallization. 

2.  Assume  that  the  formation  (crystallization)  of  a  single  stem  is  an 
elementary  process  or  elementary  reaction;  that  the  destruction 
(melting)  of  a  single  stem  is  an  elementary  process  or  elementary 
reaction;  and  that  transition  state  theory  can  be  applied  to  these 
two  elementary  processes  with  a  single  transition  state  corresponding 
to  a  relative  free  energy  maximum  or  barrier  thus  occurring  between 
each  two  integral  values  of  >/  on  a  plot  of  vs.  v. 

3.  Assume  that  only  one  stem  at  a  time  can  be  formed  or  destroyed. 

The  kinetic  problem  is  to  derive  an  expression  for  the  net  rate  Sj^(i,T) 

at  which  stems  of  length  i  (and  width  a)  pass  over  or  surmount  the  i/th  free 
energy  barrier  at  temperature  T.  The  problem  requires  consideration  of  the 


following  set  of  connected  elementary  reactions 


6 

S 


3 


6 


4  .  .  . 


where  A  is  the  rate  constant  for  the  forward  reaction  v 


1/+1  (t/  >  1)  and  B  is 


that  for  the  reverse  reaction  i/+l  -►  v(i/>l),  and  where  Aq  and  are  the 
analogous  rate  constants  for  the  i/-0  ^  j/-l  reactions.  Solution^  of  this 


problem  in  the  steady- state  approximation  gives 


S^(i,T)  - 


NoAo(A-B) 

A-B+Bj^ 


S(i.T) 


for  all  u,  where  Nq  is  the  number  of  sites  or  locations  available  for  the 
placement  of  a  first  stem.  The  total  net  rate  at  which  stems  (i.e.  the  net 
rate  including  stems  of  all  possible  values  of  i)  pass  over  the  i/th  barrier  at 
temperature  T  is  given,  for  all  i/,  by 
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STotal<T)  -  >  S(i.T) 


where  is  Che  smallest  allowed  value  of  i  which  satisfies  the  constraint 


i  > 


Af  • 


Note  that  i  is  a  discrete  variable- -the  smallest  increment  in  i  that 


can  be  made  is  the  monomer  repeat  length  To  find  Jt-^,  first  write  i  - 


2(7. 


where  m  is  a  positive  integer  and  then  realize  that,  since  i  must  exceed 
the  smallest  value  of  m  must  be  equal  to  the  smallest  integer  greater  than  ^ 

S  r  1  2a  /Af 

■J  .  Therefore,  -  |^l+E(x)Ji^,  where  x  ■  — -  and  E(x)  designates  the 


2a^/Lf- 


u 


2a. 


integer  part  of  x.  Substituting  into  the  expression  for  gives 

-  •  To  a  good  approximation,  “  T  (i.e.  x  is  sufficiently 


2a 

greater  than  1)  so  that  £-^  ~ 

CO 

Finally,  one  assiimes  that  \  S(i,T) 


4-i: 

U  ■'  i-i 


S(i,T)di;  and  the 


kinetically-determined  average  lamellar  thickness  is  then  given  by 

h  £  S(i,T)  di 

i(T)  -  - . 

S(i,T)  di 

■*1 


IV.  DETERMINATION  OF  THE  RATE  CONSTANTS 

To  obtain  expressions  for  Aq,  Bj^,  Aj^  and  B,  one  must  first  determine 
expressions  for  the  free  energy  barriers  for  the  relevant  reactions  v  ^  i/+l  (i/ 
>  0).  Let  Ej^  be  the  free  energy  barrier  to  the  destruction  of  the  first  stem; 
then  A^2^+  Ej^  is  the  barrier  to  the  formation  of  the  first  stem  in  order  that 
(A<^j^+E^)  -  Ej^  -  Let  E2  be  the  free  enrgy  barrier  to  the  formation  of 

each  subsequent  stem;  then  E  +  E2  is  the  barrier  to  the  destruction  of  each 
such  stem  in  order  that  (E+E2)  -  E2  -  E.  Now  one  does  not  know  the  free 
energy  barrier  to  the  formation  of  the  first  stem.  At  least,  one  does  know 
that  it  depends  on  what  length  i'  of  a  fully  adsorbed  stem  of  length  £ 
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actually  crystallizes  before  the  barrier  is  surmounted.  If  V  -  0,  then  none 

of  the  free  energy  of  fusion  is  released  before  the  barrier  is  surmounted,  and 

clearly,  -  2ah<7^+2h£a  and  —  abiAf.  In  general  then,  for 

0  <  £•  <2,  A.^^+Ej^  -  2aba^  +  2bio  -  abi'Af  and  Ej^  -  ab(i-i')Af.  Since  £•  is 
unknown,  a  parameter  ^  j  with  0  <  ^  <  1,  is  introduced  in  order  that  all 
possible  so-called  apportionments  of  the  free  energy  of  fusion  abiAf  between 

the  rate  constants  for  the  formation  and  destruction  of  a  first  stem  (i.e.  for 

the  forward  and  reverse  reactions  0^1)  can  be  considered.  Thus, 

A^j^  +  Ej^  -  2abcTg  +  2bia  -  ^abiAf  and  E^^  -  (1-^)  abiAf. 

Note  that  the  greater  the  amount  ^abiAf  of  the  free  energy  of  fusion  which  is 
in  fact  "apportioned"  (i.e.  the  greater  the  value  of  ^  or  i'),  the  smaller  the 
value  of  both  A^^^  +  E^^  and  E^  (for  a  given  i  and  T)  .  A  very  similar 

3 

interpretation  of  ^  has  been  discussed  recently. 

Similarly,  for  each  subsequent  stem,  let  i"  (0  <  i"  <  i)  be  the  length  of 

a  fully  adsorbed  stem  which  actually  crystallizes  before  the  barrier  to  the 

formation  of  the  stem  is  surmounted.  Then  Eo  “  2aba  -  abi"Af  and  E  +  E<,  - 

2  e  L 

ab(i-i")Af .  Define  the  apportionment  parameter  0  with  0  <  0  <  1  so  that 

E2  -  ^abOg  -  ^abiAf  and  E  +  £2  ■"  (1-0)  abiAf. 

Finally,  utilizing  transition  state  theory. 


Aq  -  h« 


-  (A^j^+Ej^+A^)/kT  -  (A^2^+E2^)/kT 


^  e 


-Ei/kT  -E,AT  -(E+E,)AT 

,A-^e‘^  ,B-^e  ^ 


where  Af"  is  the  contribution  to  each  barrier  as  a  result  of  retardations  in 


the  transport  of  a  polymer  chain  through  the  liquid  to  the  substrate  or  vice 
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B  A 

versa.  Note  that  ^  does  not  depend  on  y  and  that  does  not  depend  on  rf>  as 


required. 


V.  THE  EXTENSION  OF  THE  LAURITZEN -HOFFMAN  MODEL 

As  implied  throughout  the  above  discussion,  the  application  of  transition 

state  theory  to  the  elementary  processes  of  single  stem  formation  and 

destruction  presumes  that  there  is  a  single  relative  free  energy  maximum*^  or 

barrier  between  each  two  integral  values  of  i/  on  a  plot  of  L<t>^  vs.  v. 

Consequently,  A^j^+E^^,  Ej^,  E2,  and  E+E2  should  never  be  negative.  learly,  E^^- 

(1-^)  abiAf  and  E+E2  -  (l-^)abiAf  are  never  negative;  however,  the  expressions 

given  above  for  A^^+Ej^  and  E2  can  be  negative.  In  fact,  E2 ,  for  example,  is 

2a 

negative  for  all  i  such  that  - —  <  i  <  <»  for  a  given  Af,  (ff,  and  a  .  We 

^Af  ® 

propose  to  extend  the  LH  model  by  incorporating  into  the  model  the  assumption 
that  free  energy  barriers  cannot  be  negative,  i.e.  only  apportionments  of  the 
free  energy  of  fusion  which  result  in  a  nonnegative  barrier  will  be  allowed. 

In  order  to  incorporate  this  constraint  into  the  model,  first  note  that 
A<^j^+E^  -  2aha^  +  2bia  -  ^ablAf  is  never  negative  when  Ai^^^  is  positive  since 
then,  abiAf  <  2aha^  +  2biCT  always  holds  and  V’^biAf  <  2aha’^  +  2bia  follows. 
However,  when  is  negative,  the  expression  2aha^  +  2bia  -  ^abiAf  can  be 
negative.  The  requirement  that  A^j^  +  E^^  >  0  hold  when  A^j^  is  negative  implies 
that  one  is  not  allowed  to  apportion  all  of  the  free  energy  of  fusion  abiAf 
when  A^j^  is  negative.  If  the  amount  ^abiAf  of  the  free  energy  of  fusion  which 
is  apportioned  were  to  exceed  2aba^  +  2bia,  then  A^j^+Ej^  would  be  negative. 

The  maximum  amount  which  can  be  apportioned  is  indeed  2aha^  +  2bi<7 ,  and 
therefore  one  has,  when  A^j^  <  0, 

Ai^j^+Ej^  -  ^(2abag  +  2bicr) 

where  ^  is  an  apportionment  parameter  with  0  <  ^  <  1.  Using  " 

A<^2^  or  E^  -  (Ai^j^+Ej^)  -  A^2^  gives 

Ej^  -  £(2aba^  +  2bia)  -  (2aba^  +  2bio -abiAf)  -  abiAf- (1-^)  (2aba^+  2bia)  . 
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Observe  that  the  requirement  that  ^  0  holds  when  is  negative  is 

equivalent  to  the  physically  realistic  requirement  that  the  barrier  £2^  to  the 
destruction  of  the  first  stem  cannot  be  smaller  than  the  free  energy  increase 
(-Ai^2^)  that  occurs  upon  its  destruction.  (Note  that  abiAf -  (2abag+2bia)  - 
-A(^2^.)  Also,  this  physically  realistic  requirement  implies  that  an  adsorbed 
first  stem  cannot  completely  crystallize  before  the  barrier  to  the  formation 
of  that  stem  is  surmounted,  i.e.  that  the  upper  limit  on  i'  is  less  than  i 
when  A^2^  is  negative.  (This  upper  limit  on  i'  is  determined  later.)  For  A^^^ 

>  0,  the  expressions  A^^^  +  -  2abo^  +  2bia  -  ^abiAf  and  Ej^  -  (l-^)ab2Af 

still  hold  with  0  <  ^  <  1  and  0  <  i'  :S  i. 

At  this  point,  a  simple  change  of  variable  is  introduced  for  convenience. 
Define  \  ■  1-^  with  0  <  A  <  1. 

Now  our  approach  would  appear  to  have  introduced  another  parameter  A  in 
addition  to  if),  but  this  is  not  the  case.  To  see  this,  first  observe  that  when 
A^2^  is  positive,  the  free  energy  of  fusion  which  can  be  apportioned  is  abiAf, 
the  amount  in  fact  apportioned  is  \&abiAf,  and  the  fraction  of  the  free  energy 
of  fusion  which  can  be  apportioned  that  is  in  fact  apportioned  is  When  A^^^ 
is  negative,  the  free  energy  of  fusion  which  can  be  apportioned  is  abiAf  - 
(-A<^2^)  -  2ah(7^  +  2bi<7,  the  amount  in  fact  apportioned  is  A(2abag  +  2bia)  ,  and 
the  fraction  of  the  free  energy  of  fusion  which  can  be  apportioned  that  is  in 
fact  apportioned  is  A.  If  we  always  choose  the  same  value  for  A  and  then 

over  the  whole  range  of  values  for  A^j^,  the  fraction  of  the  free  energy  of 
fusion  which  can  be  apportioned  that  is  in  fact  apportioned  has  the  same 
value.  Let  7  denote  any  particular  value  which  is  chosen  for  both  i>  and  A, 
where  0  <  7  <  1. 

As  will  be  seen,  the  symbol  7  has  been  introduced  for  clarity. 

Recall  that  V*  ■  'J~'  expression  for  A  in  terms  of  i'  or  vice  versa  still 

needs  to  be  obtained,  and  furthermore,  equal  values  of  ^  and  A  do  not  in 
general  imply  the  same  value  of  i'.  In  our  approach,  then,  i'  depends  at 
least  on  the  sign  of  A<^2^,  but  nevertheless,  our  approach  has  only  one 
parameter,  7- -the  fraction  of  the  free  energy  of  fusion  which  can  be 
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apportioned  that  is  in  fact  apportioned- -which  is  a  constant  over  the  whole 
range  of  values  for 

In  summary,  the  barriers  in  terms  of  the  apportionment  parameter  7  are 

—  (l-7)(2abag  +  2bia) 

for  A^j^  <  0 

-  abiAf  -  7(2aba^  +  2bit7) 


A^j^  +  Ej^  —  2abag  +  2bi<7  -  7abiAf 

} 


for  A^^  >  0 


-  (l-7)abiAf 

where  we  now  observe  that  (l-7)(2ab<7g  +  2bio)  -  2aba^  +  2bia  -  7abiAf  when  Ai^^ 
-  0,  i.e.  A^^  +  E|^  is  a  continuous  function  of  i  and  Af  at  the  points  (i,Af) 
for  which  h<i)^  -  0.  Note  that  the  greater  the  value  of  the  apportionment 

parameter  7,  the  smaller  the  value  of  both  A^j^  +  E^^  and  Ej^. 

Next,  an  expression  for  t'  in  terms  of  A  is  to  be  derived.  Given  A^j^  + 

Ej^  -  (1-A)(2abag  +  2bia)  for  A^j^  <  0,  one  can  first  find  ^  when  Ai^j^  <  0  holds 

in  terms  of  A  by  equating  the  expressions 

(1-A)(2ab<7^  +  2bia)  -  2aba^  +  2bia  -  V>abiAf 


whence 


^  - 


r2a' 
_ e 

iAf 


_22l 

aAfJ 


Clearly,  equating  these  expressions  and  expressing  when  A^^^  <  0  in  terms  of 
A  is  valid  since  decreasing  Zaba^  +  2bia  by  an  amount  ^abiAf  must  be 
equivalent  to  decreasing  2aha^  +  2h£a  by  Xi2aha'^  +  2bia)  .  Note  that  the 

expression  +  aAf]  always  less  than  one  when  A^j^  is  negative.  (To  see 

this,  simply  observe  that  A^j^  <  0  implies  Zaba^  +  2bia  <  abiAf,  and  then 
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divide  both  sides  of  this  inequality  by  abiAf.)  But  ^  ^  for  all  values  of 

AijS^  so  that 

^  -221 

'^'^tiAf  ^  aAfJ  • 

Note  that  since  A  cannot  exceed  one,  the  largest  possible  value  of  i' ,  i.e. 

/2a' 

the  upper  limit  on  i'  ,  is  i  +  aAfl  ^^1  ^  mentioned  previously, 

this  upper  limit  is  indeed  less  that  i  for  A^j^  <  0. 

For  completeness,  one  can  also  find  A  when  A(Aj^  >  0  holds  in  terms  of  V  by 
equating  the  expressions 


(1-A)(2aba^  +  2bia)  -  2aba^  +  2bia  -  ^abiAf 


whence 


A  - 


[iAf  ^  aAf  J 


Clearly,  equating  these  expressions  and  expressing  A  when  A^^  >  0  in  terms  of 
Vi  is  valid  since  decreasing  2aba^  +  2bia  by  an  amount  VabiAf  must  be 

0  I 

equivalent  to  decreasing  2abag  +  2bia  by  A(2aba^  +  2bia) .  Here  again,  ^  . 

/2a'  2n  1 

Note  that  [1^  +  is  always  greater  than  one  when  A^j^  is  positive. 

In  summary,  then,  for  A^j^  <  0,  one  chooses  a  value  from  zero  to  one  for 

the  parameter  7,  whence  A  -  7,  and  then  calculates  rfi  - 

>  0,  one  chooses  a  value  from  zero  to  one  for  the  parameter  7,  whence  -  y, 

.2a;  ^  2ay-  • 

[iAf  aAfJ 


and  then  calculates  A  - 


A  -  7 


0 


for  A^j^  ^  0 
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r/>  -  y 


X 


[iAf  ^  aAf  J 


^  for  >  0 


Incidentally,  the  constraint  2aba^  +  2bia  -  ^abiAf  >  0  combined  with  0  <  <  1 

implies  that  the  inequality 

0  <  ^  <  1  the  smaller  of  1  and 
must  be  satisfied,  and  clearly  our  theory  has  satisfied  it. 

Similarly,  the  constraint  abiAf  -  A(2abtrg  +  2biCT)  >  0  combined  with  0  <  A  <  1 
implies  that  the  inequality 

0  <  A  <  the  smaller  of  1  and  — ^ - 

[iAf  ^  aAf  J 

must  be  satisfied,  and  clearly  our  theory  has  satisfied  it. 

The  approach  developed  above  can  readily  be  applied  to  incorporate  into  the 
model  the  constraint  that  E2  be  nonnegative.  Here,  E2  "  ^aba^  -  ^abiAf  can  be 
negative  when  E  is  positive,  and  E  is  always  positive  (except  when  i  -  2a^/Af, 
which  gives  E  -  0) .  The  requirement  E2  St  0  implies  that  one  is  not  allowed  to 
apportion  all  of  the  free  energy  of  fusion  abiAf.  If  the  amount  ^abiAf  which 
is  apportioned  were  to  exceed  2aba^,  then  E2  would  be  negative.  Therefore, 
one  has  E2  -  rf2aha^  where  is  an  apportionment  parameter  with  0  <  »j  <  1.  And 
E+E2  -  -2aha^  +  abiAf  +  ri2aiba^  -  abiAf  -  (l-r;)2abag.  For  convenience,  make 
the  change  of  variable  ■  1-fj  with  0  <  ^  <  1  so  that  for  all  i  and  Af 

E2  ”  (l-5)2ab<7g  and  E  +  E2  ”  abiAf  -  62aha^. 

Observe  that  the  barrier  E  +  E2  to  the  destruction  of  the  second  and  each 
subsequent  stem  cannot  be  smaller  than  the  free  energy  increase  E  that  occurs 
upon  its  destruction,  which  implies  that  an  adsorbed  second  or  subsequent  stem 
cannot  completely  crystallize  before  the  barrier  to  the  formation  of  that  stem 
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is  surmounted,  i.e.  that  the  upper  limit,  determined  below,  on  i"  is  less  than 
i. 

Given  E2  -  (1-S)2aba^,  one  can  find  ^  in  terms  of  6  by  equating  the 
expressions 

a-e)2aba^  -  2abOg  -  <^abMf 

whence 

^  fAf  • 

Clearly,  equating  these  expressions  and  expressing  ^  in  terms  of  6  is  valid 
since  decreasing  E2  with  ^  -  0  by  an  amount  ^abiAf  must  be  equivalent  to 
decreasing  it  by  S2aha^.  Note  that  the  constraint  2aba  -  ^abiAf  >  0  implies 

A 

that  the  inequality  0  <  y  <  must  be  satisfied;  since  0  <  <  1  holds,  we 

2a 

have  indeed  satisfied  this  inequality.  Also  note  that  is  always  less  than 

2a 

or  equal  to  one  since  i  >  has  been  established.  (Incidentally,  2aba^  • 

(pabiAf  >  0  does  not  imply  constraints  t  <  - —  ,  Af  <  - —  ,  or  a„  >  . ) 

^Af  e  2 

A 

Finally,  recalling  that  V  “  J~  and  substituting  above  gives  t"  -  9—^. 

In  the  special  case  -y  -  6  —  0 ,  our  model  reduces  to  the  case  ^  ^  -  0  of 

the  LH  model  which  permits  negative  barriers  for  nonzero  }/>. 

VI.  DETERMINATION  OF  THE  SIGN  OF  A^^ 

At  this  point,  one  needs  to  determine  when  A^^,  is  positive,  zero,  and 
negative.  Now  A^,  -  2aba'  +  2bfa  -  abfAf  >  0  implies  bi(2a-aAf)  >  -2aba' ;  and 

X  ®  © 

there  are  three  cases  to  consider. 

7n  -2aba' 

Case  (a):  2a  -  aAf  >  0  or  Af  <  “.  Then  the  inequality  f  >  5^2T^aAf^ 

is  always  satisfied  since  f  is  always  greater  than  zero,  and  hence  A^^^  >  0 

holds . 

Case  (b)  :  2a  -  aAf  -  0  or  Af  -  Then  A^j^  -  2abag,  which  is  always 

positive  or  zero  depending  on  a^. 
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Thus,  combining  cases  (a)  and  (b)  ,  we  have  >  0  for  all  i  when  Af  < 
(So  far,  A^j^  is  zero  only  if  both  -  0  and  Af  - 


Case  (c) ;  2a  -  aAf  <  0  or  Af  >  Then  A^,  >  0  implies  -bi(aAf-2a)  > 

Si  X 

— Thus,  when  Af  >  ^,  L<f)^  >  0  holds  for  i  <  ip,  and 
aAf 

2a 


-2ab<7'  or  i  < 
e 


Ai^j^  <  0  holds  for  f  >  iQ.  (Observe  that  as  Af  —  from  values  greater  than 

2i 

a 


iQ  -►  <».)  There  is,  however,  one  further  condition  to  consider  here. 


2a  2a 

Recall  that  i  >  has  been  established.  If  iQ  <  holds,  then  i  >  iQ 

holds  and  consequently  A^j^  <  0  would  hold  for  all  i.  To  determine  when  iQ  < 

2a' 


2a. 


_e 

Af 


2a. 


2a 


holds,  simply  write  - and  noting  that  <  1,  rearrange  this 

aAf 

inequality  to  get  Now,  if  ^  this  inequality  would  be 

e 

never  satisfied;  hence  iQ  <  never  occurs  when  a^  <  a'^.  If 
2a_ 


<  0, 

which  is 

‘"e  > 

^e'  ^0  ^ 

<  ^2 

f  1 

“  a 

k-'J' 

>  ^ 

[/%•] 

d 

occurs  when  Af  > 


Zs. 


Thus,  if  a  >  a'  and  —  <  Af 
e  e  a 


VII.  EXPRESSIONS  FOR  (T)  AND  2(T) 


If  Og  <  a^,  our  model  with  no  negative  barriers  has 


(1)  -  2abag+2bia-7abiAf 

(2)  A(^^+E2^  -  2ab<7g+2bi<7-7abiAf 

(2)  A<^^+E^  -  (l-7)(2aba;+2bia) 


for  Af  < 

a 


2a 


for  Af  >  —  and  i  <  i,, 
a  U 


for  Af  >  ^  and  i  >  i^ 
a  0 
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and  if  a  >  a'  , 
e  e 


(1)  -  2abag+2bi<7-7abiAf 

(2)  A^j^+E^  -  2abag+2bif7-7abiAf 

(2)  A^j^+Ej^  -  (l-7)(2ab<7^+2bi£T) 

(3)  A^j^+E^  -  (l-7)(2ab<7^+2bia) 


for  Af  <  ^ 
a 

for  <  Af 
a 

for  <  Af 
a 

for  Af  >  ^ 
a 


and 


and 


The  purpose  of  categories  (1),  (2),  and  (3)  will  be  seen  shortly. 

When  A^^+E^  —  2ab<7g  +  2bi<7  -  7abiAf,  E^^  —  (l-7)abiAf,  which  we  call  Case 
I. 


When  A<^2^+E^  -  (l-7)(2ab<7g  +  2bia)  ,  Ej^  -  abiAf 
call  Case  II. 

One  always  has 


E2  -  (l-tf)2aba^ 


E+E2  - 


■2aba  +  abiAf  +  E,,  -  abiAf 
e  z 


Also , 


S(i.T)  - 


A> 

1-  £  +  ^ 
^  A  ^  A 


7(2ab£7'  +  2biff)  ,  which  we 
6 


Blaho 

e 


where  f  -  e'^AT  ^ 
A  ’  A 


-(E,-E2)AT  -(A^,+E,)AT 

e  ^  ,  and  Aq  -  ^  e 


Abbreviate  c' 

abAf 
kT 


oc  2bg 
i,  "  kT 


2aba' 
_ e 

kT 


and  ^ 


2aba  2a 

— '  “  ■  a^  ’  recall  i^ 

-c  +  i.  For  Case  I, 

“^1 


2<7e 

"aT- 


wr~ 


.  2biq  -yabiAf 
kT  kT  *  kT 


c'  +  /-  (a-7)i 
■*1 
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Vfl  _  n-r)abMf  .  n  y)2  Q  e)c 

kT  kT  kT  i, 


For  Case  II, 


(I-7)  (2abag+2bi<7) 


-  (l-7)c'  +  ^  (i-7)>e 
■*1 


£..-£«  abiAf-7(2aba'+2bi<7)  (l-tf)2aba 

kT - kT  - - kT -  “  (l-a7)i-7c'-(l-«)c 


For  Case  I, 


For  Case  II, 


^NQe-*"'  (i.gC  g-ci/ip 

■  l.e<= 

n  »p\  _  _ y _ _ _ _ _ 


Sjjd.T)  -  ^  ,(1-«)C^7C'  ^-a.o,y)c!/!j 


For  any  Af  in  category  (1) ,  then. 


/^  £  Sj(i,’ 

St^Li^T)  -  r-  L  Sj(i,T)  di  and  5<1>(T)  -  -1 - 

""  1  Sj(i,T)< 


For  any  Af  in  category  (2) , 


StoLi<''>  -  ij-  {*“  *  17  j] 


- ^ 


iSj(i,T)di  +  j]  i  Sjj(£,T)d£ 


Sj(i,T)di  +  j^^Sjj(£,T)di 


For  any  Af  in  category  (3), 


n  iS,j(i,T)di 

4oLi<’'>  -  ij-  -  7=^^ - — 


|^^Sjj(i,T)di 


For  purposes  of  comparison,  the  LH  model  which  permits  negative  barriers 
has,  for  all  i  and  Af, 

-  2ab<7g  +  2bia  -  ^frabiAf  and  E2  —  2abag  -  ^abiAf 
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so  that 


and 


and 


^  -±  i 


c 


s(LH)(^_T) 


^Nq  ^-c\-c(a-^)£/£^ 

1-6*=  ^-a-^hc£/£^ 


^  n  is^^)(i.T)di 

-  r  n  s'“>(^T)dJ  and  J<“>(T)  -  -  1  - 

^  J.  S^^^(i,T)di 

1 

As  is  the  case  in  the  LH  model,  our  model  has  two  parameters.  The  most 
logical  choice  for  fl  is  5  -  7;  however,  even  with  -  7,  our  integrals  cannot 
be  evaluated  analytically.  There  seems  to  be  no  special  case  (other  than  S  - 
7  -  0)  for  which  they  could  be  evaluated  analytically.  At  this  point  then,  we 
proceed  without  setting  B  —  y. 


VIII.  EVALUATION  OF  THE  S,j,^tal^^^  2(T)--THE  VARIABLE  TRANSFORMATIONS  FOR 

THE  NUMERICAL  INTEGRATIONS 

The  required  numerical  integrations  were  easily  performed  interactively 
on  the  VAX  using  the  IMSL  subroutine  DQDAGS^  Integrals  to  be  evaluated  using 
DQDAGS  cannot  have  an  infinite  limit  of  integration.  One  way  to  proceed 
before  using  DQDAGS  is  to  make  a  change  of  integration  variable.  Although 
DQDAGS  can  integrate  functions  with  endpoints  singularities  (when  the 
endpoints  are  finite),  a  change  of  variable  which  results  in  a  transformed 
integrand  which  is  bounded  at  all  points  including  the  finite  endpoints  in  the 
new  range  of  integration,  is  preferable  to  a  change  of  variable  which  yields 
an  improper  integral  albeit  with  finite  integration  limits.  For  each  of  the 

integrals  appearing  in  S.^otal^^^’  ^Total^^^’  ^Total^^^’  ^  variable 

transformation  which  resulted  in  a  proper  integral  was  in  fact  found.  The 
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same  transformations  did  not  transform  the  corresponding  integrals  in  the 
numerators  of  and  into  proper  integrals;  however, 

the  transformed  integrands  were  of  the  form  (-In  x)f(x)  with  the  singularity 
resulting  only  from  the  factor  In  x  as  x  -►  0.  This  endpoint  singularity  could 
be  handled  by  DQDAGS . 

Consider  first  the  integral  in  The  variable  transformation 

consists  of  defining 

X  -  . 

Note  that  x(i  -*  ®)  -  0;  the  constant  i.e.  the  i  -  independent  factor, 

is  chosen  so  that  x(i  -  i^)  —  1.  Solving  for  i  in  terms 

of  X  gives  i  -  ij^Tl  -  provided  y  ^  1.  Then  di  -  - 

Furthermore,  -  e'^'  x^  ,  and 

-(l-7)ci/ii  -(1-7)0  _ 

e'  '1-e'  xso  that 


(T)-- 
^Total^^^  i 


Simplifying  gives 


(l-7)c 


r^_jci2 

Jn  u 


(J)  dx 


l-x^-T^  + 


^n)  fl 

^Total^^^  “  i  (l-7)c  Jo  _±_ 

^  ,  ..1-7 


0-1  1 

x^-yg-x^-y) 


l-x^-^  +  X 

This  is  one  of  the  integrals  that  was  evaluated  numerically  by  DQDAGS. 
Designate  the  integrand  above  as  f  j^(x)  .  Using  the  same  variable 
transformation  to  evaluate  the  numerator  of  i^^^(T)  gives 

}(l),-,v  Jo^l  [^~(i--f)c]^l^^^^  »  -^1  x)f]^(x)dx 

^  ~  fl  ~  ''  (l-7)c  fl 

Jo  f2^(x)  dx  Jo  fj^(x)  dx 

f‘^0 

Next,  using  the  same  transformation  on  the  integral  Sj(i,T)  di 
appearing  in  S^q^^j^(T)  gives 
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Sj(i.T)di  -  ^Nq 


e-c'g-(a-7)c„  , 

(l-7)c  ^ 


where 


x.  -  x(i-i!.)  -  g- (l-7)c'/(l-a) 


with  iQ  "  2ag/(l-Q)Af  as  defined  previously. 


Similarly,  the  integral  iSj(i,T)di  appearing  in  becomes 

Ji°  ^Si(i.T)  di  -  (i-,)c - ^  {^1  Jx/l^^)^^  ^  (Trfe  x)f,(x)dx} 


A  different  transformation  is  made  on  the  integral  Sjj(i,T)  di  also 


(2) 

appearing  in  Here,  define 


,(l-7)(c-c')  g-(l-7)aci/i^ 


Again  x(i  -♦  oo)  -  0;  the  constant  ^  is  chosen  so  that  x(i-iQ)  -  Xq, 

if  T 

which  is  given  above.  Solving  for  i  gives  i  -  ^  Ic-c'-  I  provided  7  ^ 

1.  Then  di  -  -  y, —  (-)  dx.  Furthermore,  g*  -g*  (^-t)  (c*c' )  ^ 

(l-7)ac  'x' 


g-ci/ii  _  ^-(c-c')/a  ^(l-7)a  ^  (l-a7)ci/ii 

Substituting  gives 


l-g-Y 

,(1-7)q 


Sjj(i,T)di- 


'  e  ' 
(l-7)ac 


.  r^o  xf 

l-e^e 

tl'] 

L<1-7)<.  1 

Jo  .  - 

[c-c' 

1 

/I  _  /  n 

l-e'"  e'^  “  J 


T -  (i)  cix 

l-gy  ^x' 
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^Nq 

r  ^0 

l-e'^e 

fc-c''] 

1 

La  J 

x(l-7)a 

(l-7)ac 

0 

l-e*^e 

c-c ' 

1  _ 1 

.-(tf-7)c  gC  g*| 

c-c' 

1  1-a-Y 

1  «  J 

1  x^^-'^>“  +  ( 

1  a  J 

1  x(l-7)a 

dx 


Designate  the  integrand  above  as  f2(x).  Thus, 


s(2)  fl  ^  ^  f«0, 

^Total''^^  "  I 


fl  >1  j-ySNo  pcQ  ^ 

—c - (i-7)«“‘  Jo 


Similarly,  the  integral  iSjj(i,T)di  appearing  in  2^^^(T)  becomes 

(l-7)ac  {  QC  -^ijo  (l-7)ac  Jo  (‘I"  x)f2(x)dx| 


Therefore , 


^5(2) 


(T)  - 


fejj°  j!/s„(i.i)di] 


ia. 

/3N, 


c(2)  /TN 

^Total^^^ 


(2) 

with  the  appropriate  expressions  for  the  integrals  and  ^jotal^^^ 
substituted  above. 

Finally,  consider  the  integral  in  Sj^^^j^CT).  The  variable  transformation 
to  be  made  on  this  integral  is 


X  ■  e 


(l-7)ac  ^-(l-7)Qci/ij^ 


Again  x(i  -*00)  -  0  and  the  constant  e^^  t)oc  chosen  so  that  x(i-ij^) 
Solving  for  i  gives  i  -  ^1-  (l-7^ac  ]  1  ’*  !•  Then  di  - 

1 _  .1.  .  r-  -ci/i, 

T—  dx.  Furthermore,  e'  "  'l-e'  "  x,e  '1- 


-  1. 


(l-7)ac  'x 
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l-g-Y 


,  and  so  that 

AN  --(1-7)(C'-K»c)p  , 

c(3)  ^  ^ ^  f  - 

^Total’'^''  i  (l-7)ac  JO 


l-g-y 


dx 


l.x(^'7)o  +  g-^c  g7(c'+gc)  j^(l-7)g 


Designate  the  integrand  above  as  f2(x).  Using  the  same  transformation  to 
evaluate  the  numerator  of  gives 

.  f^(-ln  x)fo(x)dx 

-^1  JQ  _ i _ 

Jo 


•  (3) 


(T)  -  i,  +  ., 

'  '  1  (l-7)gc 


IX.  RESULTS  AND  DISCUSSION 

A  VAX  FORTRAN  program  was  written  to  evaluate  the  required  mathematical 

expressions.  The  program  contains  the  statement  CALL  DQDAGS;  and  this  IMSL 

subroutine  performed  the  numerical  integrations .  All  calculations  were  done 

double  precision  using  the  model  parameter  values  given  in  Figure  3  of 

-8  2  2 

Reference  1;  namely,  a  -  b  -  5  x  10  cm,  a  -  10  erg/cm  ,  -  100  erg/cm  ,  T^ 

-  500  K,  Ah  -  3  X  10^  erg/ca? ,  and  Af  -  Ah(T^  ^The  average  lamellar 

thickness  calculated  from  the  IN  model  is  independent  of  ;  this  is  true  for 
our  model  only  for  Af  <  ^,  however.  Other  quantities  such  as  do 

depend  on  even  in  the  IN  model,  and  physically,^  one  expects  0  < 

In  the  case  ■  0,  our  model  is  slightly  simpler,  for  then 

A^j^  +  —  2hta  ~  yabiAf  Af  <  ^ 

-  (l-7)2bia  Ll  > 

Let  us  investigate  our  model  in  detail  for  the  case  <7^-0  first;  this  is  also 

1  2 

the  somewhat  arbitrary  choice  for  a'^  made  for  the  calculations  ’  for  the  IN 
model . 
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Given  the  parameter  values  above  and  now  with  the  choice  ^-7,  the 
calculated  average  lamellar  thickness  vs.  temperature  curves  (2  vs  T)  are 
plotted  in  Figure  1  for  the  selected  values  of  7  —0.0,  0.50,  0.75,  0.90,  and 
0.95.  Some  of  the  data  used  to  construct  these  plots  is  given  in  Table  I. 

(For  Af  <  ^,  the  average  lamellar  thickness  is  given  bj'  the  expression  for 
;j(l)(T)  given  previously  and  for  Af  >  by  the  expression  for  also 

given  previously.)  The  effect  of  7  on  2  as  a  function  of  T  is  readily 
apparent.  For  7  -  0  and  7  -  j,  2  decreases  monotonically  with  decreasing  T, 
but  for  7  “  2’  indeed  a  break  or  discontinuity  in  the  slope  cf  2  vs. 

T,  albeit  barely  discernible,  at  the  temperature  T*  for  which  Af  -  i.e.  at 
T  -  433^  K.  (This  statement  will  be  qualified  later.)  As  for  7  -  2  for  7 

-  .9,  and  .95  decreases  with  decreasing  T  for  all  T  for  which  Af  >  ^,  and 

there  is  a  break  in  the  slope  of  2  vs  T  at  T  —  T*.  Unlike  for  7  “  2’ 
plotted  2  vs  T  curves  for  7  ”  2^.  .9,  and  .95  pass  through  a  relative  minimum 

A 

at  a  temperature  for  which  Af  <  the  temperature  T^  at  which  this  relative 

3 

minimum  occurs  appears  to  increase  with  increasing  7  (for  7  “  2^.  it  occurs 

between  T  -  440  and  433-j  K  and  so  can  hardly  be  seen  on  the  plot).  Also,  over 

the  interval  T  <  T*,  2  vs  T  is  at  a  relative  maximum  at  T  -  T*;  and  2 

increases  more  rapidly  as  T  approaches  T*  both  from  values  greater  and  from 

values  less  than  T*  the  larger  the  value  of  7.  Note  that  at  least  for  all 

values  of  Af  >  — ,  2  at  a  given  T  is  larger  the  larger  the  value  of  7. 

For  comparison,  we  have  reproduced  Figure  3(b)  of  Reference  1  as  our 

Figure  2,  which  shows  the  LH  model  2  vs.  T  curves  with  ^  ^  for  the  selected 

1111 

values  of  -  0,  -j,  .9,  and  .95.  Some  of  the  data  which  we 

calculated  in  order  to  construct  these  plots  is  given  in  Table  II.  The  LH 
model  ^  -  0  curve  is  identical  to  our  7  -  0  curve.  For  Af  <  each  of  the 
LH  model  "V>  curves"  is  qualitatively  similar  but  not  quantitatively  identical 
to  its  corresponding  "7  curve"  presented  in  Figure  1.  Recall  that  the 
quantitative  difference  arises  from  the  fact  that  the  barrier  E2  has  been 
constrained  to  be  nonnegative,  i.e.  ^2  “  (i‘^)2ab<7g.  For  Af  >  ^,  however, 
the  LH  model  xp  curves  are  in  marked  contrast  to  the  7  curves;  in  particular. 


24. 


for  each  curve,  i  approaches  Infinity  asymptotically  as  Af  approaches 


2a 
V>a ' 


This  is  the  behavior  which  is  known  as  the  5i  catastrophe.  The  above  features 

of  the  LH  curves  for  <  .95  also  apply  to  LH  curves  for  .95  <  ^  <  1; 

calculations  for  the  special  case  ^  -  1  can  be  done  using  the  analytical 

expression  obtained  from  in  the  case  ^  The  curves  for  .95  <  ^  < 

1  are  similar  to  the  ^  -  .95  curve;  Figure  3(b)  of  Reference  1  gives  a  sketch 

of  the  -  1  curve,  which  exhibits  the  5i  catastrophe  at  Af  -  — . 

By  contrast,  the  7  curves  for  high  values  of  7  less  than  one  do  not 

exhibit  an  infinite  average  lamellar  thickness.  Curves  for  7  -  0.99  and  7  - 

0.998  are  presented  in  Figures  3  and  4,  respectively,  and  do  exhibit  the 

3 

features  described  previously  for  the  7  -  2^,  .9,  and  .95  curves.  Again,  for 
2^ 

Af  <  — ,  the  curves  for  7  -  .99  and  7  -  .998  are  qualitatively  very  similar  to 
LH  curves  with  V  “  -99  and  V  “  .998,  respectively.  Calculations  for  7  >  .998 
as  well  as  for  7  ■  1  apparently  cannot  be  done  using  the  expressions  for 
^(1)(T)  and  as  a  result  of  the  factor  (I-7)  appearing  in  various 

denominators . 

One  point  is  worth  emphasizing  here;  namely,  the  relationship  between  7 

0  < 

and  V-  In  both  our  model  and  the  LH  model,  i’  ^  J~<  I*nt  this  ratio  in  the  LH 
model  is  a  constant,  whereas  in  our  model 


0 


-( 


-  (ct  -  it-] 


A^..  <  0 


A(^,  >  0 


For  the  case  a'  -  0,  this  becomes 
e  ’ 


0 


-1 


2a 
^  aAf 


Af 

a 


Af 

a 
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Now,  for  any  given  say  ,  I  in  the  LH  model  is  infinite  for  all  Af  > 

J  Vja 

and  for  all  Af  >  there  is  no  finite  value  of  2  for  any  r/>  >  \f>y 

Equivalently,  a  value  of  ^  is  not  possible  for  a  chain-folded  system  for 

all  Af  >  that  is,  high  values  of  ^  do  not  lead  to  chain-folded  polymer 

crystals  at  high  enough  supercooling  according  to  the  LH  model.  Experiment, 

however,  gives  chain- folded  crystals  at  high  supercooling  with  an  average 

lamellar  thickness  that  decreases  monotonically  with  decreasing  temperature. 

As  we  have  seen,  our  one -parameter  (i.e.  y)  model  with  a'  -  0  does  reproduce 

this  high  supercooling  behavior.  And  yet,  high  values  of  }/),  i.e.  of  the  ratio 
i ' 

are  not  associated  with  our  high- supercooling  chain-folded  systems.  To 

see  this,  first  introduce  the  dimensionless  quantity  x,  where  0  <  x  <  1.  Then 

for  any  Af  -  — ,  ^  -  7  “  TX.  Since  7  cannot  exceed  one,  r/)  in  our  model 

cannot  exceed  x.  for  any  Af  >  where  x.  is  any  given  value  of  x.  But  this 

J  Xja  j 

is  exactly  what  was  found  for  rp  in  the  LH  model,  i.e.  that  a  value  of  V 

greater  than  or  equal  to  ip.  is  not  possible  for  any  Af  >  t^.  Thus,  for  Af  > 

J  Vja 

our  model,  through  the  Imposition  of  the  constraint  that  barriers  be 
nonnegative,  places  exactly  the  same  upper  limit,  on  our  ip  that  is 

predicted  for  ip  in  the  LH  model.  However,  for  Af  >  our  model,  unlike  the 

a 

LH  model,  predicts  2  vs .  T  in  qualitative  agreement  with  experiment  for  all 
values  of  our  model  parameter  7  less  than  one. 

Thus,  the  selected  calculations  done  for  our  model  indicate  that,  for  the 
case  £7^—0,  our  model  does  not  exhibit  an  infinite  average  lamellar  thickness 
for  any  value  of  7  less  than  one.  Host  importantly,  at  least  for  all  T  for 

A 

which  Af  >  our  calculations  indicate  that  our  a'  -  0  model  predicts  2  vs. 

T  curves  which  are  monotonically  decreasing  with  decreasing  T  in  agreement 

with  experiment.  That  is,  we  have  successfully  extended  the  LH  model  to 

highest  supercooling.  This  success  coupled  with  the  numerical  results  shown 

in  Figure  1  increases  our  confidence  in  using  low  values  of  7  such  as  7  -  0  as 

3 

a  first  approximation  for  convenience  in  practice.  Our  model  curves  do  show 
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an  anomaly,  i.e.  a  break-in  slope,  at  T  -  T*.  We  strongly  suspect  that  there 
is  indeed  a  break  in  slope  at  T  -  T*  because  the  relation 


2a 


Af  > 


2g 

a 


Af  < 


<  2s. 


implies  that  ^  is  discontinuous  at  Af  -  however,  we  have  not  yet 

evaluated  ^  at  Af  -  This  anomaly  is  apparently  negligible  up  to  y  values 

of  about  where  the  slope  of  2  vs.  T  has  the  same  sign  (positive)  regardless 
of  whether  the  point  Af  -  ^  is  approached  from  values  of  Af  higher  or  lower 
than  As  y  increases,  however,  this  anomaly  becomes  pronounced  with  the 

concomitant  appearance  of  a  relative  maximum  in  2  at  T  -•  T*  and  a  relative 
minimum  in  2  at  T  -  T*;  necessarily  then,  the  slope  of  2  vs .  T  as  Af 
approaches  ^  from  values  less  than  ^  becomes  negative.  We  will  refer  to 

a  o 

this  undesirable  behavior,  manifest  at  high  values  of  y,  as  the  2  anomaly. 
Unlike  the  62  catastrophe  in  the  LH  model,  the  relative  maximum  in  2  vs.  T,  as 
noted  above,  always  appears  at  Af  -  ^  for  all  values  of  y  given  that  a^  -  0. 
Despite  the  2  anomaly,  we  see  that  the  exclusion  of  negative  barriers- -the 
only  difference  between  our  model  and  the  LH  model- -has  strengthened  the 
Lauritzen-Hoffman  approach  to  polymer  cirystallization. 

One  set  of  results  with  6  f*  y  is  presented  in  Table  III.  Here,  we  see 
that  for  7  “  2  6  I,  the  calculated  2(T)  differ  only  slightly  from  the 

case  with  7  ”  2  ^  2' 

Next,  we  investigated  our  model  for  a’^  r*  0.  (Recall  that  2  for  the  LH 

model  is  independent  of  a^  and  that  our  model  is  independent  of  a^  for  Af  < 

^.)  Using  the  same  values  for  a,  b,  a,  a^,  T^,  and  Ah  as  above  and  again 

2 

with  8  -  y,  2  vs .  T  curves  for  -  0,  60,  100,  and  150  erg/cm  --each  with  7  - 

^--are  plotted  together  in  Figure  5.  Some  of  the  r*  0  data  used  to 
construct  these  plots  is  given  in  Tables  IV  and  V  (and  the  a^  -  0  data  has 
been  seen  previously  in  Table  I).  From  Figure  5,  we  see  that  2  vs.  T  is 
relatively  insensitive  to  the  value  of  for  7  ”  2’ 
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n 

erg/cm  ,  1  decreases  monotonically  with  decreasing  T,  although  for  —  100 

O 

erg/cm  ,  the  2  vs.  T  curve  is  almost  flat  near  T  -  405  K.  For  — 

150  erg/cm^,  there  is  a  relative  minimum  in  2  vs.  T  near  T  —  405  K,  and  the 
curve  passes  through  a  small  and  "diffuse"  relative  maximum  at  a  lower 
temperature .  Recall  that  one  expects  0  <  so  that  with  - 

9  9 

100  erg/cm  ,  -  150  erg/cm  may  not  be  realistic  but  is  examined  in  order  to 

explore  the  model  predictions  as  a  function  of  a'^ . 

9 

The  2  vs.  T  curves  for  a'  —  0,  60,  100,  and  150  erg/cm  --each  with  y  — 

^  --are  presented  in  Figure  6.  The  curves  pass  through  a  common  relative 
minimum  between  T  -  440  and  433^  K  (for  which  Af  <  ,  and  then  each  curve 

rises  and  passes  through  a  relative  maximvun,  that  maximum  being  relatively 
higher  and  occurring  at  higher  Af  the  larger  the  value  of  a'^.  At  each 
maximum,  there  would  appear  to  be  a  break  in  the  slope  of  2  vs.  T;  the  slope 
of  2  vs.  T  as  Af  approaches  ^  both  from  values  greater  and  from  values  less 
than  ^  is  greater  in  magnitude  the  larger  the  value  of  o'  Having  passed 
through  its  maximiun,  each  curve  decreases  monotonically  with  decreasing  T 
thereafter. 

One  should  be  careful  to  note  that  what  appears  to  be  a  break  in  the 


slope  of  2  vs.  T  when  f*  0  Ls  probably  not  a  break  in  slope; 

be  continuous  for  all  relevant  T.  Whether  a  break  in  the  slope  of  2  vs.  T 
2o 


should 


when  o'  0  as  was  presumed  true  for  a'  -  0  cannot  be 
a  e  e 


occurs  at  Af  - 

determined  conclusively  from  the  appearance  of  the  graphs,  although  the  break 
appears  to  be  absent. 

Qualitatively  similar  2  vs.  T  curves  are  obtained  for  7-0.9  and  -  0, 
60,  100  and  150  erg/cm^  as  is  shown  in  Figure  7.  Here,  the  relative  maxima 

3 

are  higher  and  "sharper"  than  the  corresponding  7  "  ^  curves,  and  they  have 
moved  to  higher  temperature.  For  y  -  0.99,  the  analogous  curves,  shown  in 
Figure  8,  exhibit  2  values  which  are  unrealistically  large  as  well  as  maxima 
which  are  extremely  "sharp". 
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The  relationship  between  7  and  ^  with  <7^  0  is  worth  emphasizing  at  this 

i ' 

point.  To  reiterate,  in  both  our  model  and  the  LH  model,  ^  ■  j-,  but  this 
ratio  in  the  LH  model  is  a  constant,  whereas  in  our  model 

f  I'  =  0 

!&(i.T)  -  \ 

[  7  A^^(i.T)  >  0 

where  the  notation  ^(i,T)  and  A^j^Ci.T)  emphasizes  here  the  dependence  of  ^  and 
on  i  and  T.  (The  T  dependence,  of  course,  enters  through  Af.)  Recalling 
the  conditions  which  govern  the  sign  of  A^j^  then  gives  when 


f!fe  ^  _2g] 

liAf  aAfJ 


\&(i.T)  - 


for  all  i  when  Af  > 


2a  r  .^e.  .] 

—  a  Lag-a^J 

^  for  i  >  iQ  when  ^  <  Hf  <  ^  (^"V] 

^  tfk  A*' 


{ 


for  i  :S  ip  when  <  Af  < 
for  all  i  when  Af  <  ^ 


and  when  a  <  a' 
e  e 


T^(i.T) 


lie  2ol.] 

iAf  ^  aAfJ 


7 


for  £  >  £q  when  Af  >  ^ 

for  i  <  when  Af  >  ^ 
u  a 

for  all  £  when  Af  <  ^ 

a 


2a' 


where  £q  - 


Af 


1- 


aAf 


Now,  on  an  i  vs.  T  curve,  one  has 
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-1 


[iAf  aAfJ 


<  0 

A^3^(i,T)  >  0 


where  the  conditions  which  govern  the  sign  of  A^j^Cl.T)  are  those  given  above 
for  A^2^(i,T)  but  with  i  replaced  by  2.  Therefore,  the  temperature  Tq  of  a 
point  (2q,  Tq)  on  an  2  vs.  T  curve  and  at  which  A^j^(2,T)  -  A^j^(2q,Tq)  -  0  is 
the  solution  to  the  following  non-linear  algebraic  equation  in  the  one  unknown 
T: 


»(2) 


(T)  -  if 


or 


A 


iSjdi+ 


Af 


f^O 
J  9. 


Sidi+ 


n 


1  - 


2g 

aAf 


If  a  >  a' 
e  e 


a 


Af  > 

a 


,  Tq  will  correspond  to  a  value  of  Af  in  the  range  —  <  Af  < 
but  if  <7g  <  a^,  Tq  will  correspond  to  a  value  of  Af  in  the  range 


Rather  than  attempt  to  solve  the  above  equation  iteratively,  one  simply 
plots  the  left-hand  side  2^^^(T)  vs.  T  and  the  right-hand  side  iQ(T)  vs.  T  on 
the  same  graph,  and  Tq  is  given  by  a  point  of  intersection  of  the  two  curves. 
Note  that  as  Af  approaches  ^  from  values  greater  than  iQ  approaches 

2^ 

infinity  and  that  decreases  monotonically  with  decreasing  T  for  Af  >  — . 

For  each  of  the  2  vs .  T  curves  with  0  presented  in  Figures  5  through  8, 

we  found  one  point  of  intersection  (2q,  Tq) ,  which  is  designated  on  each  curve 
in  the  figures  by  an  open  circle.  We  also  found  that  2^^^(T)  >  2q  holds  when 
T  <  Tq  and  that  2^^^(T)  <  2q  holds  when  T  >  Tq.  Thus,  A<^(2,T)  <  0  holds  for  T 
<  Tq  and  A«^(2,T)  >  0  holds  for  T  >  Tq.  Our  final  result  is  that,  on  an  2  vs. 

T  curve. 


Tf.  <  T  <  T° 
u  in 
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Note  that  if  the  dimensionless  quantity  x,  0  <  x  <  1,  is  again  introduced  by 

+  1 

Xj  for  some  Af  >  where  Xj  is  any  given  value  of  x.  Clearly,  one  can 

easily  calculate  i>  for  each  point  on  the  0  i  vs .  T  curves,  but  we  have 

not  yet  done  so. 

Thus ,  from  the  graphs ,  we  see  that  the  i  anomaly  becomes  more  pronounced 
but  moves  to  higher  temperature  as  y  increases  for  a  fixed  nonzero  value  of 

That  is,  although  the  relative  maximum  in  J  vs.  T  can  appear  at  some  Af  > 
when  Og  is  nonzero,  the  maximum  becomes  less  pronounced  as  it  moves  to 
lower  temperature  upon  an  increase  in  y.  Our  model,  then,  does  not  fail  at 
high  supercooling,  but  does  exhibit  anomalous  behavior  for  temperatures 
corresponding  to  values  of  Af  "Just"  greater  and  "just"  less  than  This 

undesirable  behavior  is  pronounced  for  large  values  of  y  and  is  more 
pronounced  for  larger  values  of  a'^  for  a  given  y. 

We  can  easily  rationalize  mathematically  how  our  calculated  2  vs .  T 
curves  can  rise  with  decreasing  T  for  some  Af  >  ^  when  is  nonzero.  Recall 
that  the  expression  for  l^^^(T),  namely 


j  so  that,  unlike  the  case  —  0,  ^  can  exceed 


writing  Af 


_  22. 

xa’ 


if>  -  7X 


raa' 


[la 


/  °fSj(i,T)df  +  r  iS,j(i,T)di 

r^o  r 

Sj(i,T)di  +  Sjj(i,T)di 


contains  two  different  integrands  Sj(i,T)  and  Sjj(i,T).  Depending  on  a^,  y, 
and  T,  the  contribution  of  the  integrals  involving  Sj(i,T)  to  may 

outweigh  the  contribution  of  the  integrals  involving  Sjj(i,T),  and  in  some 
cases,  our  calculations  show  that  to  a  very  good  approximation 

fS.(f.T)di 

l^  '  (T)  ~  -  with  £q  approaching  infinity. 

J  °  3.(i,T)di 
1 
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But  this  is  our  expression  for  for  the  interval  Af  i  and  the 

results  of  our  calculations  using  have  been  found  to  differ  little 

from  results  using  i.e.  the  LH  theory.  Not  unexpectedly  then, 

can  increase  with  decreasing  T  for  some  Af  >  We  note  that  the 

numerator  of  Sj(i,T),  like  the  numerator  of  S^^\i,T),  contains  the  factor  Aq 
-  e  ^  e  form  of  which  has  been  associated  with^  increases 

in  J  with  decreasing  T. 


X.  CONCLUSIONS 

Thus,  the  1  anomaly  is  apparently  connected  to  the  expression  A^^  "*■  ^1  “ 

2aba^  +  2bio  -  yabiAf  even  when  the  maximvun  in  2  vs .  T  occurs  at  a  temperature 

for  which  Af  exceeds  Our  results  with  -  0  clearly  Indicate  that  the  2 

anomaly- -and  in  part  the  62  catastrophe  of  the  LH  theory- -are  associated  with 

the  interval  Af  <  ^  and  are  thus  connected  to  the  expression  A^j^+Ej^  -  2aha^  + 

2b£a  -  yabiAf.  Even  when  high  values  of  7  or  0  are  considered  unrealistic  as 

3 

has  been  elucidated  recently,  however,  there  is  no  guarantee  that  the  LH 
theory  as  well  as  our  extension  of  it  has  not  failed  to  incorporate  an  as  yet 
unknown  constraint  or  feature  which  would  improve  the  model  results  at  high  7 
values.  For  example,  high  7  values  may  be  unrealistic,  but  the  2  values  for 
high  7  from  an  improved  model  may  simply  be  unrealistically  large  but 
nevertheless  monotonically  decreasing  with  decreasing  T  for  all  T.  Further 
work  along  this  line  would  probably  be  mostly  of  theoretical  interest  rather 
than  essential  for  use  in  practice. 

Although  the  2  anomaly  remains  in  our  model,  we  have  shown  that  the  62 
catastrophe  of  the  LH  theory  is  in  part  related  to  the  failure  to  exclude 
negative  barriers.  Furthermore,  our  model  is  successful,  for  it  shows  that 
the  Lauritzen -Hoffman  approach  to  polymer  crystallization,  subject  to  the 
exclusion  of  negative  barriers,  is  consistent  with  experimental  behavior  at 
very  high  supercooling.  We  conclude  that  we  can  safely  extend  our  version  of 
the  LH  theory  to  treat  systems  interacting  with  an  external  electric  field. 
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FIGURE  CAPTIONS 

Figure  1.  Plots  of  Average  Lamellar  Thickness  vs.  Temperature  for  7  -  0, 
0.90,  and  0.95,  each  with  “  0.  See  text  for  a,  b,  a, 
a^,  T^,  and  Ah  which  are  the  same  for  Figures  1  through  8. 

Figures  1  through  8  have  $  —  y. 

At  T  -  433J  K  (i.e.  Af  -  ^),  -  0. 

For  T  >  433^  K,  A^j^  i  0  and  1^-7. 

For  T  <  433^  K,  :S  0  and  ^  -  7(^). 

Figure  2.  Plots  of  Average  Lamellar  Thickness  vs.  Temperature  for  ^  -  0, 
113 

3’  2’  4’  0.95  reproduced  from  the  Lauritzen- Hoffman 

Model  (Reference  1);  plots  are  independent  of  <7^. 

Figure  3.  Plot  of  Average  Lamellar  Thickness  vs.  Temperature  for  7  -  0.99 
and  —  0.  Again,  for  T  >  433^  K,  —  7  and  for  T  <  433^  K, 

^  “  T'fSf)- 

Figure  4.  Plot  of  Average  Lamellar  Thickness  vs.  Temperature  for  7  -  0.998 
and  o'  -  0.  Again,  for  T  ^  433^  K,  ^  -  7  and  for  T  <  433j  K, 

(if) 

Figure  5.  Plots  of  Average  Lamellar  Thickness  vs.  Temperature  for  “  0, 

2  1 

60,  100,  and  150  ergs/cm  ,  each  with  7  “  2*  Each  open  circle 
designates  the  point  (ig,  Tg)  at  which  A^j^(i,T)  -  0.  For  T  >  Tg, 

A^i  >  0  and  ^  —  7.  For  T  i  Tn.  <  0  and  ^  —  7  f — —  +  . 

^  u  i 

Figure  6.  Plots  of  Average  Lamellar  Thickness  vs.  Temperature  for 

2  3 

-  0,  60,  100,  and  150  ergs/cm  ,  each  with  7  “  4-  Again,  each 
open  circle  identifies  the  temperature  Tq  (see  Figure  5). 

Plots  of  Average  Lamellar  Thickness  vs.  Temperature  for 

2 

a’  -  0,  60,  100,  and  150  ergs/cm  ,  each  with  7  -  0.90.  Again, 


Figure  7 . 
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Figure  8 . 


each  open  circle  identifies  the  temperature  Tq  (see  Figure  5). 

Plots  of  Average  Lamellar  Thickness  vs.  Temperature  for 
<7g  -  0,  60,  100,  and  150  ergs/cm^,  each  with  7  -  0.99. 

For  -  0,  60,  100,  and  150  ergs/cm^,  Tq  -  433  ^  K,  432.2  K, 
432.1  K,  and  432.0  K,  respectively  (see  Figure  5). 
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TABLE  CAPTIONS 

Table  I.  Average  Lamellar  Thickness  as  a  Function  of  Temperature  for 

7  “  0.  2’  0.90,  each  with  —  0  and  0  —  y.  See  Figure  1. 

Table  II.  Average  Lamellar  Thickness  as  a  Function  of  Temperature  for 
^  and  0.90  reproduced  from  the  Lauritzen-Hoffman  Model 
(Reference  1),  each  with  ^  and  independent  of  a^.  See 
Figure  2 . 

Table  III.  Average  Lamellar  Thickness  as  a  Function  of  Temperature  for 
7  -  ^  -  1,  and  <7^-0.  See  text  for  the  usual  values  of 

a ,  b ,  <7 ,  <7^  ,  T  “ ,  and  Ah . 

’  *  *  e '  DU 

Table  IV.  Average  Lamellar  Thickness  as  a  Function  of  Temperature  for 

2  1 

t7g  -  60,  100,  and  150  ergs/cm  ,  each  with  1  ~  2-  See  Figure  5. 

Table  V.  Average  Lamellar  Thickness  as  a  Function  of  Temperature  for 
a'^  -  60,  100,  and  150  ergs/cm^,  each  with  7  -  0.90.  See 
Figure  7 . 


Average  Lamellar  Thickness  (A) 


Average  Lamellar  Thickness  (A) 


Average  Lamellar  Thickness  (A) 


Avarage  Lamellar  Thickness  (A) 


Temperature  (K) 


Pigurs  5 


Average  Lamellar  Thickness  (A) 


Table  I.  Average  Lamellar  Thickness  (X)  vs.  Temperatiire  (K) 


TEMP.  (K)  psi=Gamma=0  Gamma=1/2 


/ 

TEMP.  (K)  Gamma=0.90 


485.000 

234.383 

235.303 

485.000 

236.013 

480.000 

1  78.390 

1 79.781 

480.000 

180.939 

475.000 

1 44.660 

146.556 

475.000 

148.300 

470.000 

122.074 

124.507 

470.000 

127.027 

465.000 

105.867 

108.867 

465.000 

1 12.435 

460.000 

93.652 

97.253 

460.000 

102.279 

455.000 

84.105 

88.342 

455.000 

95.475 

450.000 

76.429 

81.344 

450.000 

91.672 

445.000 

70.1 15 

75.762 

445.000 

91.275 

440.000 

64.826 

71.267 

440.000 

96.119 

435.000 

60.328 

67.641 

435.000 

1 12.616 

430.000 

56.451 

63.528 

430.000 

117.625 

425.000 

53.072 

59.481 

425.000 

109.730 

420.000 

50.100 

55.988 

420.000 

103.882 

415.000 

47.463 

52.941 

415.000 

99.316 

410.000 

45.105 

50.259 

410.000 

95.563 

405.000 

42.984 

47.877 

405.000 

92.353 

400.000 

41,064 

45.744 

400.000 

89.529 

395.000 

39.316 

43.821 

395.000 

86.992 

390.000 

37.718 

42.077 

390.000 

84.676. 

385.000 

36.251 

40.484 

385.000 

82.538 

380.000 

34.897 

39.023 

380.000 

80.545 

375.000 

33.644 

37.676 

375.000 

78.673 

370.000 

32.480 

36.429 

370.000 

76.905 

365.000 

31.396 

35.270 

365.000 

75.225 

360.000 

30.382 

34.188 

360.000 

73.622 

355.000 

29.433 

33.176 

355.000 

72.087 

350.000 

28.540 

32.225 

350.000 

70.612 

345.000 

27.700 

31 .329 

345.000 

69.190 

340.000 

26.907 

30.484 

340.000 

67.817 

335.000 

26.157 

29.683 

335.000 

66.486 

330.000 

25.446 

28.924 

330.000 

65.194 

325.000 

24,772 

28.201 

325.000 

63.938 

320.000 

24.130 

27.513 

320.000 

62.714 

315.000 

23.518 

26.855 

315.000 

61.519 

310.000 

22,935 

26.226 

310.000 

60.352 

305.000 

22.377 

25.624 

305.000 

59.210 

300.000 

21.843 

25.045 

300.000 

58.090 

295.000 

21.332 

24.489 

295.000 

56.992 

290.000 

20.841 

23.953 

290.000 

55.913 

285.000 

20.369 

23.437 

285.000 

54.853 

280.000 

19.915 

22.938 

280.000 

53.809 

275.000 

19.479 

22.456 

275.000 

52.782 

270.000 

19.057 

21.990 

270.000 

51.769 

265.000 

18.651 

21.537 

265.000 

50.770 

260.000 

18.258 

21.099 

260.000 

49.784 

255.000 

17.878 

20.673 

255.000 

48.810 

250.000 

1  7.51 1 

20.259 

250.000 

47.847 

245.000 

17.155 

19.856 

245.000 

46.895 

240.000 

1  6.809 

19.463 

240.000 

45.953 

235.000 

16.475 

19.081 

235.000 

45.021 

Table  II.  Average  Lamellar  Thickness  (X)  vs.  Temperature  (K) 


TEMP  (K)  t-H  Psi=1/2  LH  Psia0.90 


485,000 

235.785 

237.166 

480.000 

180.224 

182.177 

475.000 

146.926 

149.552 

470.000 

124.780 

128.225 

465.000 

109.027 

1 13.507 

460.000 

97.290 

103.129 

455.000 

88.251 

95.962 

450.000 

81.124 

91.560 

445.000 

75.412 

90.139 

440.000 

70.789 

93.098 

435.000 

67.037 

105.777 

430.000 

64.009 

160.924 

425.000 

61.610 

Oo 

420.000 

59.786 

415.000 

58.519 

410.000 

57.832 

405.000 

57.800 

400.000 

58.577 

395.000 

60.458 

390.000 

64.019 

385.000 

70.494 

380.000 

82.999 

375.000 

1 12.171 

370.000 

232.547 

000 

oo 

360.000 

355.000 

350.000 

345.000 

340.000 

335.000 

330.000 

325.000 

320.000 

315.000 

310.000 

305.000 

300.000 

295.000 

290.000 

285.000 

280.000 

275.000 

270.000 

265.000 

260.000 

255.000 

250.000 

245.000 

240.000 

235.000 


Table  III 


Average  Lamellar  Thickness  (^)  vs.  Temperatvire  (K) 


TEMP  (K) 

Theta=1 

-  495.000 

675.848 

490.000 

485.000 

230.877 

480.000 

475.000 

142.184 

470.000 

465.000 

104.542 

460.000 

455.000 

84.037 

450.000 

445.000 

71.460 

440.000 

435.000 

63.333 

430.000 

425.000 

56.368 

420.000 

415.000 

50.779 

410.000 

405.000 

46.332 

400.000 

395.000 

42.690 

390.000 

385.000 

39.639 

380.000 

375.000 

37.03B 

370.000 

365.000 

34.779 

360.000 

355.000 

32.796 

350.000 

345.000 

31.035 

340.000 

335.000 

29.454 

330.000 

325.000 

28.022 

320.000 

315.000 

26.716 

310.000 

305.000 

25.516 

300.000 

295.000 

24.405 

290  000 

285.000 

23.373 

280.000 

275.000 

22.407 

270.000 

265.000 

21.501 

260.000 

255.000 

20.646 

250.000 

245.000 

19.836 

240.000 

235.000 

19.067 

Table  IV 


Average  Lamellar  Thickness  (S)  vs.  Temperature  (K) 


TEMP.  (K) 


0.5  //  60  0.5  //  100  0.5  //  150 


485.000 

235.303 

235.303 

235.303 

480.000 

179.781 

179.781 

179.781 

475.000 

146.556 

146.556 

1  46.556 

470.000 

124.507 

124.507 

124.507 

465.000 

108.867 

108.867 

108.867 

460.000 

97.253 

97.253 

97.253 

455.000 

88.342 

88.342 

88.342 

450.000 

81.344 

81.344 

81.344 

445.000 

75.762 

75.762 

75.762 

440.000 

71.267 

71.267 

71.267 

435.000 

67.641 

67.641 

67.641 

430.000 

64.735 

64.735 

64.735 

425.000 

62.454 

62.454 

62.454 

420.000 

60.723 

60.743 

60.743 

415.000 

59.214 

59.577 

59.564 

410.000 

57.306 

58.874 

59.005 

405.000 

54.856 

58.337 

58.984 

400.000 

52.149 

57.582 

59.533 

395.000 

49.469 

56.411 

60.296 

390.000 

46.971 

54.852 

60.919 

385.000 

44.708 

53.035 

61.120 

380.000 

42.683 

51.095 

60.800 

375.000 

40.874 

49.136 

60.003 

370.000 

39.252 

47.220 

58.842 

365.000 

37.791 

45.385 

57.434 

360.000 

36.466 

43.648 

55.882 

355.000 

35.253 

42.015 

54.261 

350.000 

34.136 

40.486 

52.625 

345.000 

33.100 

39.056 

51.009 

340.000 

32.131 

37.719 

49.436 

335.000 

31.219 

36.468 

47.918 

330.000 

30.358 

35.295 

46.462 

325.000 

29.541 

34.194 

45.072 

320.000 

28.766 

33.159 

43.746 

315.000 

28.028 

32.183 

42.485 

310.000 

27.324 

31.262 

41.286 

305.000 

26.652 

30.390 

40.145 

300.000 

26.009 

29.564 

39.059 

295.000 

25.392 

28.778 

38.025 

290.000 

24.799 

28.031 

37.040 

285.000 

24.229 

27.318 

36.101 

280.000 

23.681 

26.637 

35.204 

275.000 

23.152 

25.985 

34.346 

270.000 

22.641 

25.360 

33.526 

265.000 

22.147 

24.760 

32.740 

260.000 

21.669 

24.184 

31 .986 

255.000 

21.206 

23.628 

31.263 

250.000 

20.756 

23.093 

30.567 

245.000 

20.320 

22.576 

29.898 

240.000 

19.896 

22.076 

29.253 

235.000 

19.484 

21.593 

28.631 

Table  V.  Average  Lamellar  Thickness  (X)  vs.  Temperatiore  (K) 


TEMP.  (K)  0.90  //  60  0.90  //  100  0.90  //  150 


485.000 

236.013 

236.013 

236.013 

480.000 

180.939 

180.939 

180.939 

475.000 

148.300 

148.300 

148.300 

470.000 

127.027 

127.027 

127.027 

465.000 

1 12.435 

1 12.435 

1 12.435 

460.000 

102.279 

102.279 

102.279 

455.000 

95.475 

95.475 

95.475 

450.000 

91.672 

91.672 

91.672 

445.000 

91.275 

91.275 

91.275 

440.000 

96.1 19 

96.119 

96.119 

435.000 

1 12.616 

112.616 

112.616 

430.000 

174.477 

176.464 

176.575 

425.000 

201.468 

297.598 

429.985 

420.000 

167.419 

248.786 

364.401 

415.000 

144.1 10 

205.816 

292.407 

410.000 

128.691 

177.388 

245.551 

405.000 

1 17.747 

157.689 

213.621 

400.000 

109.503 

143.241 

190.563 

395.000 

103.012 

132.137 

173.106 

390.000 

97.724 

123.283 

159.389 

385.000 

93.301 

1 16.014 

148.284 

380.000 

89.520 

109.905 

139.076 

375.000 

86.231 

104.670 

131.286 

370.000 

83.324 

100.111 

124.586 

365.000 

80.723 

96.086 

118.740 

360.000 

78.366 

92.490 

113.578 

355.000 

76.210 

89.245 

108.971 

350.000 

74.218 

86.289 

104.820 

345.000 
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83.577 
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70.621 
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97.602 

335.000 

68.974 

78.739 

94.427 
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67.408 

76.559 

91.485 

325.000 

65.912 

74.510 

88.746 
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72.574 
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70.740 

83.773 
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68.993 
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67.325 

79.347 

300.000 
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65.728 
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55.678 
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53.448 
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54.763 

63.695 

255.000 
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53.558 

62.244 
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48.233 

52.380 
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240.000 
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